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Abstract 



The large transverse momentum muon pair and W boson productions 
in the quark initiated coherent diffractive processes at hadron colliders are 
discussed under the framework of the two-gluon exchange parametrization 
of the Pomeron model. In this approach, the production cross sections are 
related to the small- x off-diagonal gluon distribution and the large- x quark 
distribution in the proton (antiproton). By approximating the off-diagonal 
gluon distribution by the usual gluon distribution function, we estimate the 
production rates of these processes at the Fermilab Tevatron. 
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I. INTRODUCTION 



In recent years, there has been a renaissance of interest in diffractive scattering. These 
diffractive processes are described by the Regge theory in terms of the Pomeron (IP) exchange 
The Pomeron carries quantum numbers of the vacuum, so it is a colorless entity in QCD 
language, which may lead to the "rapidity gap" events in experiments. However, the nature 
of the Pomeron and its interaction with hadrons remain a mystery. For a long time it had 
been understood that the dynamics of the "soft Pomeron" was deeply tied to confinement. 
However, it has been realized now that much can be learned about QCD from the wide 
variety of small-x and hard diffractive processes, which are now under study theoretically 
and experimentally. 

On the other hand, as we know that there exist nonfactorization effects in the hard 
diffractive processes at hadron colliders 0-0. First, there is the so-called spectator effect 
||, which can change the probability of the diffractive hadron emerging from collisions intact. 
Practically, a suppression factor (or survive factor) " Sp" is used to describe this effect ||. 
Obviously, this suppression factor can not be calculated in perturbative QCD, and is now 
viewed as a nonperturbative parameter. Typically, the suppression factor Sf is determined 
to be about 0.1 at the energy scale of the Fermilab Tevatron 0. Another nonfactorization 
effect discussed in literature is the coherent diffractive processes at hadron colliders 0, in 
which the whole Pomeron is induced in hard scattering. It is proved in |J that the existence 
of the leading twist coherent diffractive process is associated with a breakdown of the QCD 
factorization theorem. 

Recently, we have shown that these coherent diffractive processes can be studied in the 
two-gluon exchange model and have calculated the diffractive J/ip and charm jet 



production rates at hadron colliders |10|,|TT|. The two-gluon exchange parametrization of the 
Pomeron were previously used to describe the diffractive processes in photoproduction at ep 
colliders An important feature of this model recently demonstrated is the sensitivity 

of the production cross section to the off-diagonal parton distribution functions in the proton 
0. 



As sketched in Fig.l, the color-singlet two-gluon system (representing the Pomeron) is 
emitted from one hadron and interacts with another hadron by scattering off its partons 
in a hard scattering process. This process is a coherent diffractive process. The processes 



calculated in [|T0| , p~Tj1 are the gluon induced processes, i.e., the parton (emitted from the 
upper hadron) involved in the hard process is the gluon. Here, in this paper, the processes 
that will be calculated are the quark initiated processes. Following the method introduced 
in |10|,11|, we calculate the massive muon pair and W boson production in the coherent 
diffractive processes at hadron colliders in the leading logarithmic approximation (LLA) of 
QCD. As shown in Fig.l, the massive muon pair is produced via a time-like virtual photon 
7*. So, we only need to calculate the process of a virtual photon production, pp — > q^*P + X, 
and the virtuality of the photon M 2 is equal to the invariant mass of produced muon pair. 

Because of the quark participation in these processes, we expect that the production cross 
sections of these processes are sensitive to the large- a; quark distribution in the proton. As 
has been shown in [pTO , lTj| , the gluon initiated coherent processes are sensitive to the small- a; 
gluon distribution in the proton. However, all of these coherent processes are sensitive to the 
off-diagonal gluon distribution in the proton, which is the typical property of the two-gluon 
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exchange parametrization of the Pomeron model calculated previously and here. 

The diffractive production of heavy quark jet at hadron colliders has also been studied by 
using the two-gluon exchange model in Ref. ||12|| . However, their calculation method is very 
different from ours []. In their calculations, they separated their diagrams into two parts, 
and called one part the coherent diffractive contribution to the heavy quark production. 
However, this separation can not guarantee the gauge invariance fpTlfl . In our approach, 
we follow the definition of Ref. |1, i.e., we call the process in which the whole Pomeron 
participants in the hard scattering process as the coherent diffractive process. Under this 
definition, all of the diagrams plotted in Fig. 2 for the partonic processes qp — > q^*p and 
qp — > q'Wp contribute to the coherent diffractive production. 

In addition, we must emphasize that in this paper we only consider the diffractive muon 
pair (or W) production in the "Pomeron Fragmentation region", i.e., M| ~ 4A;|, where M\ 
is the invariance mass of the diffractive final states and fcy is the transverse momentum of 
the muon pair. 

The rest of the paper is organized as follows. In Sec. II, we give the cross section formulas 
for the partonic processes qp — > 7*gp and qp — > W ± qp in LLA QCD. We use the Feynman 
rule method |14[ in the calculations. The numerical results are given in Sec. Ill for diffractive 
massive muon pair and W boson productions at large transverse momentum at the Fermilab 
Tevatron. The conclusion is given in Sec. IV. 



II. LLA FORMULA FOR THE PARTONIC PROCESS 

The Feynman diagrams of the partonic processes are the same for the 7* production 
and the W boson production except the difference of the electroweak vertex in these two 
processes. So, in the following, we will mainly focus on the calculation of the 7* production 
process. The W boson production cross section will be obtained by the similar method. 

As sketched in Fig.l, the cross section for the diffractive muon pair production at hadron 
colliders (pp at the Tevatron) can be formulated as 

da{pp — > n + n~p + X) = J dxida(qp — > /i + ^qp)f q (x 1 , Q 2 ), (1) 

where x± is the longitudinal momentum fraction of the antiproton carried by the incident 
quark of flavor q. f q {xi,Q 2 ) is the quark distribution function in the antiproton, and Q 2 is 
the scale of the hard process. da{qp — > fi + ^~qp) is the cross section for the partonic process 
gp -> n + n-qp. 

In the diffractive partonic process qp — > 7*gp, at the leading order of perturbative QCD, 
there are four diagrams shown in Fig. 2. The two-gluon system coupled to the diffractive 
proton is in a color-singlet state, which characterizes the diffractive processes in perturbative 
QCD. Due to the positive signature of these diagrams (color- singlet exchange), the real part 
of the amplitude cancels out in the leading logarithmic approximation. To evaluate the 
imaginary part of the amplitude, we must calculate the discontinuity represented by the 
crosses in each diagram of Fig. 2. 



1 For detailed discussions and comments, please see [11] 



The diagrams of Fig. 2 show that the process qp — > *y*qp is by crossing related to the 
photoproduction process 7*p — > qqp in ep collisions. But the virtualities of the photons in 
these two processes are not the same. In the photoproduction process (DIS) the virtual 
photon is space-like, while in the massive muon production process it is time-like. This 
relation between the above two processes is similar to the relation between the Drell-Yan 
process at hadron colliders and the deep inelastic scattering process at ep colliders. 

In the following calculations, we express the formulas in terms of the Sudakov variables. 
That is, every four-momenta ki are decomposed as 

ki = a { q + (3ip + k iT , (2) 

where q and p are the momenta of the incident quark and the proton, q 2 = 0, p 2 = 0, and 
2p ■ q = W 2 = s. Here s is the cm. energy of the quark-proton system, i.e., the invariant 
mass of the partonic process qp — > ^*qp. o.i and A are the momentum fractions of q and p 
respectively, k^ is the transverse momentum, which satisfies 

k iT -q = 0, k iT -p = 0. (3) 

All of the Sudakov variables for every momentum are determined by using the on-shell 
conditions of the momenta of the external-line and crossed-line particles in the diagram. 

The Sudakov variables associated with the momentum k are determined by the on-shell 
conditions of the outgoing 7* momentum k + q and the light quark momentum u — k. So, 
the variables a k and j3 k satisfy the following equations: 



a k M 2 - k 2 = a k (l + a k )M 2 x , 

A = Ml^Mk, (4) 



where fcy is the transverse momentum of the virtual photon 7*, M 2 is the virtuality of the 
photon 7* (the invariant mass squared of the muon pair), and M\ is the invariant mass 
squared of the diffractive final system (including the muon pair and the light quark jet), i.e. 

M\ = (q + uf. (5) 

Because we derive the differential cross section da/dt at t = 0, in the calculations the 
momentum transfer squared of the diffractive process qp — > ^*qp is set to be zero, i.e., 
u 2 = t = 0. This identity, together with the above equation Eq. @, gives the Sudakov 
variables for momentum u as, 

M 2 

a u = 0, f3 u = -A u 2 T = 0. (6) 
Similarly, for the momentum /, the Sudakov variables are 

Jjrp 

OLi = - — , 

2(kr,l T )-(l-p k )l* 
fa = , for Diag. 1, 2, 

a k s 

M 2 - I 2 

A = -, for Diag. 3, 4, (7) 
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where {kx, It) is the 2-dimensional product of the vector k T and It- We can see that f3i is 
not the same for these four diagrams, as it is in the case of the diffractive charm jet and 



J/ip production processes flT0| , |TT 



Among these variables, there are two small parameters, 



(3 U « 1, « 1, (8) 

which are the basic expansion parameters in the following calculations. 

With the variables introduced above, the differential cross section formula for the partonic 
process qp — > ^*qp can be written as 

da. dM\dk\da k . k\ — a k M 2 ^ j- ^: 



dt lt=0 16s 2 1Qtt 5 MV^ k ^ ' ~ K> ' M 



5(a k (l + a k ) + — "* )J2\A\\ (9) 



x lv± x 

where A is the amplitude of the process qp —>■ ^qp. We know that the real part of the 
amplitude A is zero, and the imaginary part of the amplitude A(gp — > 7*gp) for each 
diagram of Fig. 2 has the following general form, 

ImA = C F I jJt^F x u{u - k)T^u(q), (10) 
J (l T ) 

where Cp = 2/9 is the color factor for each diagram. T M represents some 7 matrices including 
one propagator. F in the integral is defined as 

F=^-g 2 s ee g f(x',x";l T ), (11) 

where 

(12) 

where the function G(x', x"\ k T ) is the so-called off-diagonal gluon distribution function ||T3 



Here, x' and x" are the momentum fractions of the proton carried by the two gluons. It is 
expected that for small x, there is no big difference between the off-diagonal and the usual 
diagonal gluon densities [T5|] . So, in the following calculations, we estimate the production 
rate by approximating the off-diagonal gluon density by the usual diagonal gluon density, 
G(x',x";Q 2 ) w xg(x,Q 2 ), where x = xp = M x /s. 

In the amplitude Eq. (0), we see that the leading logarithmic contribution comes from 
the terms in T M which are proportional to l\. So, we can expand Y ^ in terms of l T , and take 
the leading order terms (l T ), and neglect the higher order terms. 

Futhermore, we note that in the integral of Eq. (|1(]) the l T terms in V ^ coming from 
all diagrams must be canceled out by each other. Otherwise, their net sum (order of l T ) 
will lead to a linear singularity as l T — > when we perform the integration over l T . The 
linear singularity is not proper in QCD calculations. So, we first examine the behavior of 
the amplitude at the order of l T , i.e., in the limit of l T 0. In this limit, the Y ^ for each 
diagram has the following form, 
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p(l) = _p(2) 



r (3) _ _ r (4) _ jf_ (l o\ 

To obtain the above result, we have used the following equations, 

u{u - k){i- ff) = 0, 

M<l) = 0- (14) 

From Eq. (0), we can see that the contributions of at the order of l T from the four 
diagrams are canceled out by each other. There is no contribution to the amplitude Eq. (|10D 
at this order. This is what we expected as mentioned above. 

From the above analysis, we see that only after summing up all of the four diagrams 
contributions it can give a gauge invariant amplitude for the partonic process qp — > q^p. 
So, the separation of the diagrams according to Ref. [0 is not correct. 

At the next order expansion of T^, l T , the evaluation is much more complicated, but 
nevertheless straightforward. We first give the expansion result for the propagators. To the 
order of l T , they are 

1 1 

J_ l + a k -(3 k )l 2 T 4(k T ,l T ) 2 _ 2(k T ,l T ) 
92 M\ { + a k M\ + {a k M x f a k M 2 x J ' 

9s = —92- (15) 
a k 

Apart from the propagator expansion, the 7 matrices in also contain the l\ terms. 
They mostly come from the Sudakov variables (expressed in l\) of momentum Z, i.e., ai, 
(3i, and the slasher fr- Furthermore, the 2-dimensional product (fey, It) also contributes l\ 
terms. After integrating over the azimuth angle of It, we get the following result, 

J d%(k T ,l T ) 2 = I [ dl 2 T k T l 2 Tj 
J d%(k T , l T ) h=\J dl 2 T AtI t - (16) 

To simplify our calculations, the contributions from the four diagrams are added together, 
and decomposed into several terms as follows. 

The terms (defined as a-term) coming from the slasher typ in the 7 matrices, for all these 
four diagrams to being summed up together, are 

^ - -4Hr * (17) 

The terms (6-term) from o.\ = —jtt are 
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r?> = Jr[(i + a*W, - (is) 



The terms (c-term) from the propagator expansion are 



Ti ^-v^ {l+a * )s2( ^^ + «m h >>- (19) 

The last terms come from those which are proportional to "(A/r, It) j/r" ■ They are (d-term) 

= / A lt 2 1+ 2 afc K^ * y (20) 

\ivi x ) a k 

Adding up all of the above a, b, c, d terms, we get the amplitude squared for the 
diffractive process qp — > ^*qp as, after averaging over the spin and color degrees of freedoms, 

ER 2 = 162Q ' Q ae ^\t5 6 (— ) 2 [2(1 + <* k )M 2 M x - a k (M x y - 2(M 2 ) 2 ](xg(x, Q 2 )) 2 . 
9 { M x) a k 

(21) 

To obtain the above result, we have only taken the leading order contribution, and neglected 
the higher order contribution which is proportional to f3 u = M\j s. The factorization scale in 
the gluon density is very important because we know that the parton distributions at small 
x change rapidly with this scale. In Ref. ||12| , they used different scales for their two parts of 
the diagrams for the partonic process. However, as discussed in the above calculations, all 
of the four diagrams must be summed together to give a gauge invariant amplitude for the 
coherent diffractive production at hadron colliders. So, the scales of the four diagrams must 
be the same, for which we choose it to be Q 2 . The separation in |12 is not gauge invariant. 

Finally, we get the differential cross section for the partonic process qp —>■ ^*qp in the 
LLA of QCD, 

|t =0 = dM x dk T da k — -[xg(x, Q )\ (- 



M 2 [2((l + a k )M 2 M 2 x - a k (M 2 x ) 2 - 2(M 2 ) 



■[5{a k - ax) + S(a k - a 2 )\, (22) 
(M 2 X )^(M 2 X - M 2 ) 2 - Ak\M 2 x 

where ot\^ are the solutions of the following equations, 

, kl-aM 2 

a(l + a) + ^2 = °- ( 23 ) 

With this differential cross section formula for the partonic process qp — > ^qp, we can 
get the differential cross section for the partonic process of massive muon pair production 
qp ^ + ^qp, 

dcr(qp — > H + fi~qp) , a da(qp — > j*qp) . 

dM 2 dt It= ° = 3ttM 2 It '* = °' 

For the W boson production processes, the differential cross section formula da/dt can 



(24) 



be easily obtained from Eq. ( |22|) by making the following replacements, 



M -> M w , ae 2 -> 2^, (25) 



2 



^ ~4tt 
where g w = G F M^/y/2. 



7 



III. NUMERICAL RESULTS 



With the cross section formulas given in last section, we can calculate the cross section 
of the diffractive production at the hadron level. However, as mentioned above, there exist 
nonfactorization effects caused by the spectator interactions in the hard diffractive processes 
in hadron collisions. Here, we use a suppression factor Ts to describe this nonfactorization 
effects in the hard diffractive processes at hadron colliders H . At the Tevatron, the value of 
Ts may be as small as Ts ~ 0.1 jH| . That is to say, the total cross section of the diffractive 
processes at the Tevatron may be reduced down by an order of magnitude due to the 
nonfactorization effects. In the following numerical calculations, we adopt this suppression 
factor value to evaluate the diffractive production rates of massive muon pair and W boson 
at the Fermilab Tevatron. 

In the numerical calculations, we take the input parameters as follows, 

a = 1/128, M w = 80.33 GeV, G F = 1.15 x 1(T 5 GeV- 2 . (26) 

The scales for the running coupling constant and the parton distribution function are set 
to be the same. For the numerical calculations, the scale is set to be Q 2 = M 2 + k\ for 
massive muon pair production (here we still use a = 1/128 for this process approximately) 
and Q 2 = for the W boson production. For the parton distribution functions, we select 



the GRV NLO set 16 



The numerical results for the massive muon pair production at large transverse momen- 
tum in the coherent diffractive processes at the Fermilab Tevatron are shown in Fig. 3 and 
Fig. 4. In Fig. 3, we plot the double differential cross section cPa^ / dM 2 dt\ t= o as a func- 
tion of the lower bound of the transverse momentum of the muon pair /cTmin, where we 
set M 2 = 5 GeV 2 . In Fig. 4, we plot the double differential cross section as a function of 
the lower bound of the momentum fraction x\, Ximin- This figure shows that the dominant 
contribution to the cross section comes from the region of x% ~ 10 _1 . Comparing this result 
with that of the diffractive J/ip and charm jet production cross sections, we find that the 
dominant contribution region of x\ to the cross section of muon pair production here is some 
orders of magnitude larger than that of J/ip and charm jet production. This is understand- 
able, because the cross sections for the diffractive processes at hadron colliders depend on 
the parton distribution functions as the following form, 

daoc/^Q 2 )^^ 2 )) 2 , (27) 

where X\ is the longitudinal momentum fraction of the proton (or antiproton) carried by 
the incident parton (quark or gluon). In the diffractive J/ip and charm jet production 
processes, the incident parton is the gluon, so the cross section is sensitive to the small- x 
gluon distribution function in the proton. In the diffractive massive muon pair and W boson 
production processes calculated here, the incident parton is the quark, so the cross section 
is sensitive to the large- X quark distribution in the proton. However, the dependence of the 
cross sections on the second factor (g(x,Q 2 ) 2 is the same for these two types of processes 
(the quark induced and the gluon induced processes). This feature is due to the two-gluon 
exchange model calculation for the diffractive processes. 

The numerical results on the diffractive W boson production at the Fermilab Tevatron are 
plotted in Fig. 5 and Fig. 6. In these two figures, we only plot the diffractive W + production 
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cross section in the process pp — > W + pX. W~ production in (pp — > W~pX) process and 
W ± production in (pp — > W^pX) can be obtained similarly. In Fig.5, we plot the differential 
cross section da w / dt\ t =o as a function of the lower bound of the transverse momentum k Tmin . 
In Fig. 6, we plot the dependence of the cross section on Xi min , where we set k Tmin = 5 GeV. 
From this figure, we see that the dominant contribution comes from the region of X\ > 10 _1 . 
This is because the large mass of the W boson requires the value of x\ larger than that for 
the low mass muon pair production. 

For the possibility about the experimental measurement of the coherent diffractive pro- 
duction, we think that some other mechanisms may be important in perturbative QCD 
beyond the two-gluon exchange model, such as quark pair exchange. We note that quark 
pair exchange is suppressed for the case of valence quarks by the factor (3 U = M\j s and 
corresponds to the secondary Reggeons (p, u, f, . . .) exchange. So, the 'valence' and 'sea' 
quark contributions may be different for the quark pair exchange processes. Work along 
this way is in preparation. Also, the survival factor and the off-diagonal parton distribu- 
tion function may cause uncertainties as well. So, we will leave the discussions about the 
experimental measurement in a forthcoming paper. 



IV. CONCLUSIONS 

In this paper, we have derived the formula for the quark initiated coherent diffrac- 
tive processes at hadron colliders in perturbative QCD by using the two-gluon exchange 
parametrization of the Pomeron model. We have shown that production cross sections 
of massive muon pair and W boson are related to the off-diagonal gluon distribution and 
large-a; quark distribution in the proton (antiproton). To estimate the production rates for 
these processes, we have approximated the off-diagonal gluon distribution by the usual gluon 
distribution function in the proton. 

By now, we have studied the coherent diffractive processes at hadron colliders by using 



the two-gluon exchange model including the J/ijj production [10], the charm jet production 



rjfl , and the massive muon pair and W boson productions. All of these processes have a 
common feature: there is a large energy scale associated with these processes to guarantee 
the application of perturbative QCD, i.e., for J/ip production, m c for the charm jet 
production, M 2 for the massive muon pair production, and Mw for the W boson production. 
And another common feature is the sensitivity to the off-diagonal gluon distribution function 



in all these processes (for detailed discussion see Wjl)- Therefore, these diffractive production 
processes will open a useful window for the study of the off-diagonal gluon distribution 
function in the proton. 
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Figure Captions 



Fig.l. Sketch diagram for the diffractive virtual photon 7* and W boson productions at 
hadron colliders in perturbative QCD. 

Fig. 2. The lowest order perturbative QCD diagrams for partonic process qp — > ^*qp. The 
crosses in each diagram represent the discontinuity we calculated in the evaluation of the 
imaginary part of the amplitude. 

Fig. 3. The double differential cross section d 2 a^ /dM 2 dt\ t=0 at M 2 = 5 GeV 2 for the massive 
muon pair production at the Fermilab Tevatron as a function of /crmin- ^Tmin is the lower 
bound of the transverse momentum of the muon pair. 

Fig. 4. The double differential cross section d 2 a m / 'dM 2 dt\t=o at M 2 = 5 GeV 2 as a function 
of xi m i n , where we set A>rmm = 5 GeV. xi min is the lower bound of the momentum fraction 
of the proton carried by the incident quark. 

Fig. 5. The differential cross section da w /dt\ t=0 for the W + boson production in the diffrac- 
tive process pp — > W + pX at the Fermilab Tevatron as a function of A/rmin- 
Fig. 6. The differential cross section da w / dt\ t =o as a function of rci min , where we set fo^mm — 
5 GeV. 
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